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We present model atmospheres and theoretical spectra for X-ray bursters. Our models include the effects of 
Compton scattering on free electrons. The atmospheres have compositions that are mixtures of hydrogen, helium, 
and iron. For our models the ratio of color temperature, T c , to effective temperature, T e //, is in the range 1.2 - 
1.8. This ratio depends on T e //, surface gravity, and iron abundance. For fixed T e // and surface gravity, models 
with non-zero iron abundance exhibit values of T c /T e ff that are slightly lower than pure hydrogen-helium models. 



1. Introduction 

Discovery of the first X-ray burster was re- 
ported by Grindlay et al. (1976) [1] and by Belian 
et al. (1976) [2]. We believe that these objects are 
low mass binary stars which contain the accret- 
ing neutron star. Since we do not observe X-ray 
pulsations, the magnetic field of the neutron star 
is lower than 10 10 Gs. Weak magnetic field of 
such intensity allows for the existence of accre- 
tion disc. Basic features of X-ray burst are: fast 
rise ~ 1 s and decay time ~ 3 - 100 s, luminos- 
ity near the peak about 10 38 erg/s and energy per 
burst 10 39 ergs (cf. review article by Joss & Rap- 
paport 1984[3]). The idea that energy source of 
a burst is the thermonuclear flash was proposed 
by Woosley and Taam (1976) [4] and Maraschi & 
Cavaliere (1977) [5]. This model quite well ex- 
plains many global features of X-ray bursters. 

Initially neutron star spectra were assumed to 
be blackbody spectra, but this assumption is not 
valid when the atmosphere was dominated by 
scattering on free electrons. In such case spectra 
of outgoing photons are shifted towards higher en- 
ergy and the maximum is shifted in the same di- 
rection both for the coherent Thomson, and non- 
coherent Compton scattering. In electron scat- 
tering stellar atmosphere, photons are created 
by thermal processes below the photosphere, in 



so called thermalization layer. In this layer lo- 
cal temperature T is much higher than the effec- 
tive temperature T e ff, therefore maximum of the 
Planck function B V (T) occurs in higher energy 
than the maximum of B u (T e ff). Consequently, 
the blackbody temperature Tbb as obtained from 
fitting of continuum spectra, can be much larger 
than the true T e f f . Temperature Tb b is usually 
referred to as the color temperature T c . 

In principle, the ratio T c /T e ff can reach ar- 
bitrary large values in Thomson scattering stel- 
lar atmospheres, provided that true absorption in 
that atmospheres decreases to zero (Madej 1974 
[6]). Compton scattering in stellar atmospheres 
reduces T c /T e ff by unknown factor, but still this 
ratio remains larger than unity. Model atmo- 
sphere computations with account of Compton 
scattering, and exact T c /T e ff determination, are 
very difficult. 

Ebisuzaki, Hanawa & Sugimoto (1984) [7] esti- 
mated T c /T e ff w 1.4 by use the simple opacity 
formulae. Model atmospheres computations by 
London, Taam & Howard (1986) [8] showed, that 
this ratio was in the range from 1.2 to 1.7 and 
depended on few parameters: T e //, surface grav- 
ity log g and iron abundance. London, Taam & 
Howard (1986) [8] concluded that the solar iron 
abundance does not influence significantly the ra- 
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tio T c /T e ff. On the other hand, Ebisuzaki & 
Nomoto (1986) [9] claimed that absorption by iron 
decrease this ratio. 

2. Model equations 

Our model equations describe transfer of ra- 
diation subject to constraints of hydrostatic and 
radiative equilibrium. We assume the equation of 
state of ideal gas in local thcrmodynamical equi- 
librium (LTE). The equation of transfer includes 
f-f absorption from all ions, b-f absorption from 
iron ions, and Compton scattering terms. The 
equation of transfer is as follows: 
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where: B v is the Planck function, J v the mean 
intensity, e„ is the ratio of true absorption to co- 
efficient of total opacity, r„ the monochromatic 
optical depth, and \x denotes cosine of zenithal 
angle. 

Details of the above equation of transfer, 
and Compton scattering redistribution functions 
<&i(i/,z/) and $2{i / ,i /l ), and linearization proce- 
dure, are described in earlier papers (see Madej 
1991 [10]; Madej & Rozanska 2000a[ll], b[12]; 
Madej, Joss & Rozanska 2003[13]). 

3. Computations 

We calculated our models with the following as- 
sumptions: equation of state in LTE, hydrostatic 
and radiative equilibrium, plane-parallel geome- 
try, non rotating star. Our equations allow for 
large photon-electron energy exchange at the time 
of Compton scattering. Details of the computer 
code and linearization scheme were extensively 
described in earlier papers quoted in previous sec- 
tion. 

We calculated grid of models each of them com- 
puted on the 144 standard optical depth points 
(12 points per the decade) from 10 4 to 10 -8 and 
901 wavelength points. We computed few mod- 
els of hot neutron star atmospheres with effective 
temperatures T e ff—1, 2, 3xl0 7 K, and surface 
gravities g from the critical gravity up to 10 15 in 
cgs units. Models assume H and He in solar num- 
ber abundance Nh c /Nh = 0.11 and iron with 



number abundance Nf £ /Nh = 3.7 x 10 5 (ap- 
proximately solar value) and 100 times greater. 

4. Results 

We claim that our computations can measure 
both the ratio T c /T e ff and its slight changes in- 
fluenced by small differences of T e ff, logg and 
iron abundance. We obtain models of burster at- 
mospheres where the constancy of the bolometric 
flux is almost everywhere better than 1%, and 
only in few optical depth points approches 2%. 

Tables 1 and 2 show values of T c /T e ff com- 
puted over wide range of parameters T e ff, logg 
and iron abundance. We can see that for mod- 
els with pure H-He atmospheres, these ratios are 
higher than for atmospheres with nonzero iron 
abundance, of course for the same T e ff and log g. 
For higher metal abundance these ratios decrease. 
This is because increase of metal abundance di- 
minish the importance of Compton scattering, 
and T c /T e ff ratio is than determined by the re- 
distribution of outgoing flux by nongray metal 
absorption. 



Table 1 

T c /T e ff ratios for different parameters. 
Left part of the tables present models with 
N He /N H = 0.11 and N Fe /N H = 10~ 5 , and right 
part N He /N H = 0.11, and N Fe /N H = 10" 3 



T eff 


logs 


Tc/T eff 


T eff 


logs 


Tc/T eff 




14.8 


1.77 




14.8 


1.25 


3 x 10 7 


14.9 


1.56 


3 x 10 7 


14.9 


1.24 




15.0 


1.44 




15.0 


1.24 




14.1 


1.67 




14.0 


1.41 


2 x 10 7 


14.5 


1.41 


2 x 10 7 


14.5 


1.38 




15.0 


1.27 




15.0 


1.27 




12.9 


1.51 




12.8 


1.27 


1 x 10 7 


14.5 


1.22 


1 x 10 7 


14.5 


1.30 




15.0 


1.25 




15.0 


1.36 



What is interesting, for the lowest T e ff = 10 7 K 
these ratios exhibit local mimimum in some inter- 
mediate logg. Tendencies, which we observe are 
due to iron opacity and the Compton scattering. 
For the lowest surface gravity Compton scatter- 
ing dominates over thermal absorption, therefore 
shift of the peak spectra is the highest. In the 
presence of iron, thermal absorption is greater 
than in case of H-Hc atmospheres. Consequently, 
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Table 2 

Table shows ratios T c /T e ff for pure H-He model 
atmospheres (Madej, Joss & Rozahska 2003) 



T eff 


log g 






14.8 


1.79 


3 x 10 7 


14.9 


1.69 




15.0 


1.61 




14.1 


1.75 


2 x 10 7 


14.5 


1.43 




15.0 


1.29 




12.9 


1.75 


1 x 10 7 


14.5 


1.29 




15.0 


1.33 



the T c /T e ff is lower for atmospheres containing 
iron. 

In Figure 1 we show the outgoing photon spec- 
tra for various parameters. We compare black- 
body spectra with the spectra of H-He atmo- 
sphere and with the atmosphere containing iron 
in Figure la. We can see that the blackbody spec- 
tra (solid line) have the peak in lower energy and 
the BB spectrum is softer than the spectrum of 
the model X-ray burster. In lower photon en- 
ergies we obtain excess of the outgoing flux as 
compared to the blackbody spectra both for H-Hc 
models and atmospheres containing some amount 
of iron. Figures lc and Id present spectra of at- 
mospheres with iron abundance 100 times greater 
than the solar abundance, and with the solar iron 
abundance, respectively. For the latter iron abun- 
dance we note, that the iron ionisation edge ap- 
pear in emission for two lowest T e ff. Profiles of 
iron absorption lines (fundamental scries of H- 
like, and He-like ions) can be very complicated 
(see Fig. lb). In some models we obtain inverse 
emission in iron line cores, which strictly corre- 
spond to temperature inversion reproduced in our 
model atmospheres with Compton scattering. 

We stress here, that we did not obtain values 
of T c /T e ff larger than 2.0 in any of our model 
atmospheres of X-ray burst sources. 
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Figure 1. Various models of neutron star atmospheres. 

a) solid line - blackbody spectrum with T = 10 7 K; dotted line - pure H-He log g — 15, T e ff = 10 7 K; 
dashed line N Fe /N H = 3.7 x 1CT 5 , log 5 = 15, T eff = 10 7 K; 

b) resonant line profile of helium-like iron; 

c) N Fe /N H = 3.7 x 1CT 3 dotted line logg = 15, T eff = 10 7 K; solid logg = 15, T eff = 2 x 10 7 K; dashed 
logs = 15, T eff = 3 x 10 7 K; 

d) N Fe /N H = 3.7 x 1CT 5 dotted line log^ = 15, T eff = 10 7 K; solid log 5 = 15, T eff = 2 x 10 7 K; dashed 
log 5 = 15, T e ff = 3 x 10 7 K. 



